Role of volcanism and impact heating in mass extinction climate shifts

This study investigates the mechanisms underlying the varied climate changes witnessed during mass extinctions in the Phanerozoic Eon. Climate shifts during mass extinctions have manifested as either predominant global cooling or predominant warming, yet the causes behind these occurrences remain unclear. We emphasize the significance of sedimentary rock temperature in comprehending these climate shifts. Our research reveals that low-temperature heating of sulfide leads to global cooling through the release of sulfur dioxide (SO2), while intermediate-temperature heating of hydrocarbons and carbonates releases substantial carbon dioxide (CO2), contributing to global warming. High-temperature heating additionally generates SO2 from sulfate, further contributing to global cooling. Different degrees of contact heating of the host rock can lead to different dominant volatile gas emissions, crucially driving either warming or cooling. Moreover, medium to high-temperature shock-heating resulting from asteroid impacts produces soot from hydrocarbons, also contributing to global cooling. Large-scale volcanic activity and asteroid impacts are both events that heat rocks, emitting the same gases and particles, causing climate changes. The findings elucidate the critical role of heating temperature and heating time in understanding major climate changes during mass extinctions.


Effects of heating temperatures on gas production
My analyses indicate that limestone samples with high organic carbon content exhibited significant CO 2 emissions within the temperature range of 700-800 °C under slow heating at a rate of 10 °C/min (Fig. 2, represented by a red bar in Fig. 4).Additionally, SO 2 gas was released from sulfate sources, namely gypsum and anhydrite, at temperatures ranging approximately from 1000 to 1300 °C (Fig. 2, depicted as a blue bar in Fig. 4).Whereas, SO 2 gas was released from sulfide within the temperature range of 400-800 °C17 (a blue bar in Fig. 4).Combining these findings with previously published data from 1-, 10-, 100-, and 1000-min heating experiments (illustrated by red and blue dots in Fig. 4) 11 , these experimental data are categorized into three groups based on the sources of their formation: SO 2 formed from sulfide under low-temperature heating, CO 2 from hydrocarbons and carbonates under intermediate-temperature heating, and SO 2 from sulfate under high-temperature heating (1000-1300 °C in 1-5-min heating, 550-650 °C in 100-year heating) (Fig. 4).The distribution of data in each group is represented by the Arrhenius equation, demonstrating the relationship between production temperature and heating time (refer to "Methods").This alignment with the Arrhenius equation indicates that the experimental data conform to the chemical reaction role, as confirmed by the fit to the Arrhenius equation.
These experimental processes and conditions marked by slow heating and low oxygen are suitable for simulating volcanism rather than asteroid impacts.Consequently, the findings suggest that substantial amounts of SO 2  4 .Open histograms in H-A depict species extinction percentages for 2060-2080 CE 53 .The abbreviations O, F, D, G-L, P, T, J, K-Pg, and H-A correspond to the following time periods: Ordovician, Frasnian, Devonian, Guadalupian-Lopingian transition, Permian, Triassic, Jurassic, Cretaceous-Paleogene boundary, Holocene-Anthropocene.These extinction percentage data are directly comparable due to the use of similar methods, including the conventional method and substage intervals 4 .The numbers 1 to 5 denote the five major mass extinctions.Each silhouette represents a representative vertebrate animal from the respective age.The end-Jurassic mass extinction event is associated with the Morokweng impact crater formation in South Africa 3 .SST: presence of SST data.
and CO 2 gases are generated during the volcanic contact metamorphism of sedimentary rocks in Large Igneous Provinces (LIPs) through both low and high-temperature heating.Additionally, medium-temperature heating induces a significant release of CO 2 .
The heating temperatures employed in the brief laboratory experiments (ranging from a few seconds to a few days) were translated into equivalent temperatures aligned with the 100-year heating time scale, pertinent to sill contact metamorphism and derived from the Arrhenius equation (refer to the Supplementary Information in Kaiho et al. 11 ).The maximum heating time of 100 years is based on a numerical model depicting the thermal evolution of a 100-m thick sill 18 .We postulate that during combustion, decomposition reactions involving calcium carbonate, pyrite (sulfide), and sulfate lead to the formation of carbon dioxide and oxidized sulfur species, respectively, depending on the heating temperature.These reactions were incorporated into the kinetics  gas from samples containing 1 wt% gypsum (sulfate) and 1 wt% artificial anhydrite (sulfate).The curve labeled 44 amu represents the temperature for CO 2 release, while the curve labeled 64 amu (b) corresponds to the temperature for SO 2 release from gypsum and anhydrite.The molar ratio of SO 2 /g gypsum or anhydrite to CO 2 /g limestone is 4%, which corresponds to a weight ratio of 6%.

Effects of heating temperatures on soot production
In the current analyses, soot formation was predominantly observed between 800 and 1000 °C, with each maximum temperature maintained for 1 min, in limestone and mudstone powder samples under a slow heating rate of 8 °C/min, at a low air pressure of 2 × 10 -4 Pa (Fig. 3a,c).Raman analysis confirmed the presence of soot on the inside wall of the ampoules for the black materials (Fig. 3d).However, prolonged heating resulted in the disappearance of soot, with each maximum temperature maintained for 10-1000 min (Fig. 3b).
This study demonstrates that to initiate soot formation, a temperature range of 800-1000 °C is required under slow short-term heating and closed low air pressure conditions.Krestinin et al. 19 , who conducted experiments on the isothermal pyrolysis of acetylene and benzene with a pyrolysis time of 0.17 s, observed an increasing trend of soot formation from 900 to 1100 °C.Moreover, Yoshihara et al. 20 showed that polyaromatic hydrocarbons (PAHs) and soot were formed at temperatures of 1500-1800 °C and 1800-1900 °C, respectively, within a time frame of 10 −3 to 10 −2 s under shock-heated conditions.He et al. 21demonstrated that abundant soot was formed at temperatures between 1500-1900 °C and at 2 and 4 atm shock temperature and pressure during 1.5 × 10 −3 s.The model calculation of Chen et al. 22 indicates soot formation under gas temperatures of 1100-1400 °C during 10 −5 s.
To summarize, significant soot formation necessitates short-term heating, lasting less than a few minutes.The required temperature range is between 800-1000 °C under low air pressure conditions, resembling the process of volcanic eruptions.Similarly, temperatures between 900-1400 °C under atmospheric air pressure conditions during 10 −5 to 10 −1 s, and 1500-1900 °C under shock-heated conditions during 10 −3 s resulting from rapid heating like an asteroid impact, also facilitate soot formation.These soot data distributions in the production temperature/heating time figure also fit the Arrhenius equation (Fig. 4).It is noteworthy that the asteroid impact scenario allows for soot formation across a wider temperature range.Contact metamorphism is the long-term heating, www.nature.com/scientificreports/which cannot form soot as shown (Fig. 4).Volcanism can produce soot through the short-term heating process of volcanic activity.Consequently, the impact of volcanic soot on climate changes is relatively minor (Fig. 4).

Heating temperatures and heating times for SO 2 , CO 2 , and soot production
Figure 4 illustrates the relationships between heating temperatures and heating times for the production of SO 2 , CO 2 , and soot.The data points in the figure conform to the distribution pattern predicted by the Arrhenius equation.This consistent pattern suggests that lower temperatures and longer heating times can produce SO 2 , CO 2 , and soot similarly to higher temperatures and shorter heating times, aligning with the principles of the Arrhenius equation.For instance, soot production at 1800-1900 °C during 10 −3 s corresponds to that at 1000 °C during 1 min (Fig. 4).The variations in heating conditions are as follows: slow heating under low oxygen conditions corresponds to contact metamorphism by magma intrusion like a sill, whereas rapid heating under atmospheric oxygen conditions corresponds to impact heating.Figure 4 represents the primary outcome of this research.Leveraging the insights from Fig. 4 and incorporating the coronene index, elaborated upon in the next section, this paper has generated Figs. 5 and 6, which automatically present the conclusions.High-temperature heating resulting from the asteroid impact primarily   11 and the present study), SO 2 from sulfate (from the present study), and soot from hydrocarbon (based on Krestinin et al. 19 [isothermal pyrolysis of acetylene, benzene, and diacetylene], Yoshihara et al. 20 under a shock-heated condition, He et al. 21under 2 and 4 atm shock temperature and pressure, model calculation of Chen et al. 22 , and the present study).The data points are represented by red and right-side blue bars, and black dots for 1-min heating, all obtained from the present study.The term "Low oxygen" corresponds to a pressure of 2 × 10 -4 Pa.The colored areas on both sides of the frame, which encompass the data plots, represent the fitting of the experimental data to the Arrhenius equation.The Arrhenius equation elucidates the relationship between production temperature and heating time.The presence of experimental data points within these colored areas indicates how well they conform to the predicted behavior based on the Arrhenius equation.
The production temperature on the horizontal axis indicates the heating temperature in each experiment.
Vol:.( 1234567890) www.nature.com/scientificreports/released SO 3 and minor amounts of SO 2 (3% or more in the total amount of SO 2 and SO 3 ) from sulfates, leading to limited global cooling [23][24][25][26] .However, soot originating from hydrocarbons caused global cooling and played a role in the K-Pg mass extinction 15,27 .
The anomalies in the coronene index and sea surface temperature (SST) can be attributed to volcanic or impact events (Table 1).While the temperature at the impact point exceeds 10,000 °C, the heating temperature decreases to a few hundred degrees moving from the center of the crater outward.The middle temperature scale represents the temperature during 5 min of heating based on Kaiho et al. 11 and experimental data from this paper.The bottom temperature scale represents the temperature heated by sill according to Aarnes et al. 18 , lasting for 100 years.
Intermediate coronene index values (0.2-0.8), indicating medium-temperature heating (700-1000 °C for 1 min), were observed in volcanic-generated mass extinctions at the Guadalupian-Lopingian transition in South  4 .In H-A, eight dots are from Kaiho 53 , with the six open dots representing projections for the near future (2060-2080 CE).Soot amount and surface temperature anomaly at the end-Jurassic were calculated 53 and obtained using Figure 5 of Kaiho and Oshima 27 .Species extinction percentages for 2060-2080 CE are estimated based on 16 cases of global warming, environmental pollution, and forest devastation as causes 53 .The production temperature scale aligns with that in Fig. 5.
High coronene index values (0.8-1.0) were found at the end of the Frasnian in France and Belgium (including GSSP; average 0.90 [SD 0.04]) 10 , as well as the impact-generated mass extinction at the K-Pg boundary in Haiti and Spain (average 0.88 [SD 0.03]) 15 (Fig. 5).
For the four events except for G-L, coronene index values are distinctly concentrated to low at End-T, medium at end-P, and high at end-F and K-Pg, resulting in no overlapping between the low, medium, and high values.While G-L has a wide range from 0.4 to 0.9, it does not plot in the low-temperature range.
The global surface temperature anomalies during each mass extinction were summarized in Kaiho 4 .Based on the surface temperature anomalies and the coronene index, these five mass extinctions can be grouped into three categories: (i) global cooling and low coronene index observed at the end of the Triassic, (ii) global warming and intermediate coronene index observed at the G-L transition and the end of the Permian, and (iii) global cooling and high coronene index at the end of the Frasnian and K-Pg boundary (Fig. 5, Table 1).

Determining factors heating temperature
The presence of extremely high temperatures, surpassing 10,000 °C, at the center of impact craters explains the observed high coronene index at the K-Pg boundary 15 .In contrast, the heightened heating temperature at the end-Frasnian can be attributed to plume magma volcanism within Large Igneous Provinces (LIPs), deep sills, and older-age magma (Table 2).These volcanic events, characterized by explosive carbonatite-alkaline magmatism and associated kimberlite-type eruptions, exhibit temperatures around 1400 °C in the Late Devonian, as determined from the age/temperature curve of Herzberg et al. 30 .Covering extensive areas with estimated diameters ranging in the thousands of kilometers, such intense volcanic activity has the potential to generate a high fraction of coronene, a seven-ringed PAH, through the significant heating of large volumes of sedimentary hydrocarbons present in older strata 10 .
However, low-temperature heating at the end-Triassic mass extinction is attributed to a shallower mantle source, shallower sill depths, and younger mantles, exhibiting temperatures around 1370 °C in the end-Triassic (Table 3).Medium-temperature heating at the end-Permian mass extinction is due to a deep mantle source, medium sill depths, and around 1370 °C mantle temperature.
The presence of high coronene index values cannot be solely attributed to wildfires.Previous studies have indicated that the abundance of coronene and benzo[ghi]perylene, which are derived from petroleum sediments, can be associated with hydrothermal volcanism off the coast of California 29 .Thus, volcanism can contribute to the formation of six to seven-ring PAHs, including coronene.Although forest fires frequently occur during the deposition of sediments in seas, it is important to note that during non-event periods, the coronene index values are typically less than 0.2 [9][10][11][12]15 . Hig-temperature wildfires, such as the recent ones witnessed in Australia with temperatures reaching 1200 °C, can contribute to the formation of coronene.Sedimentary rocks found in seas contain PAHs derived from numerous wildfires occurring in various regions dominated by forests and grasslands.The temperature of a typical grassland wildfire is around 300 °C, leading to a low coronene index when PAHs from different types of wildfires mix.The ratio of the sum of benzo[e]pyrene, benzo[ghi]perylene, and coronene to phenanthrene (5-to 7-ring PAH index) is used as an indicator of heating events [9][10][11]31 .The 5-to 7-ring PAH index is significantly elevated in heating events and mass extinctions compared to their background values (non-extinction periods) [9][10][11] .While the coronene index values during the end-Triassic mass extinction are similar to the background values, the 5to 7-ring PAH index is higher.This indicates that the average heating temperature during the mass extinctions, including the end-Triassic low-temperature volcanic event, was higher than the average temperature of wildfires.

The impact of heating temperature on global climate
The correlation between coronene index values and gas production temperature has been analyzed using the heating experimental data on PAHs from Norinaga et al. 28 (Fig. 5).This correlation suggests that the grouping of mass extinctions based on coronene index values and global surface temperature anomalies can be associated with heating temperatures/time for gas and soot production.
Global climate changes appear to be influenced by the temperature, volume, and composition of sedimentary rocks in Large Igneous Provinces (LIPs), providing a mechanism to explain the diverse climate changes resulting from LIP volcanism 11,34 .Soot formation is observed during short-term heating, characteristic of impacts, while long-term heating in volcanic settings typically leads to the production of SO 2 and CO 2 unless there is a sudden release of heated hydrocarbons through eruptions.These SO 2 and soot must reach the stratosphere to form sulfuric acid and soot aerosols inducing global cooling, a process facilitated by the high energy of volcanic eruptions and asteroid/comet impacts.When soot formed from impact target rocks was introduced into the stratospheric SO 2 from the Deccan volcanism, the soot blocked sunlight, preventing the formation of sulfuric acid from SO 2 due to reduced sunlight.
Large Igneous Provinces (LIPs), with diameters ranging from 1000 to 2000 km, exhibit similarities in rock types and proportions, primarily composed of carbonates, evaporites, and clastic rocks 1,7,8,[34][35][36][37][38][39] (Table 3).The composition of the impact target rocks at Chicxulub also includes carbonates, evaporites, and clastic rocks 40 .These suggest that there is likely minimal bias in how emitted gases, soot, and coronene are recorded between different events.The gases and soot released from LIPs and impact target rocks exist in the amount necessary to cause global climate changes 4 , as shown in Svensen et al. 34 for the Siberian LIP, Kaiho et al. 11 for CAMP, and Kaiho and Oshima 27 for the Chicxulub impact.Many sills are present, as shown in Table 3 (thickness, see note) for all Lips.
The heating temperature significantly influences the amount of released SO 2 and CO 2 , resulting in either cooling or warming effects.The sulfate content in evaporites may vary; nevertheless, sulfate does not dissolve at low to medium temperatures, impeding the release of SO 2 gas from sulfate.In cases such as the K-Pg and end-F, the high heating temperature generated SO 2 gas from sulfate rocks abundant in sulfate.
Sulfides are present in substantial amounts in all sedimentary rocks, approximately 100 times that of sulfate, resulting in the emission of SO 2 gas at low, medium and high temperatures across various LIPs.The common occurrence of hydrocarbons and carbonates in target rocks for impacts and LIPs leads to significant CO 2 generation at both moderate and high temperatures 11 .Consequently, low-temperature heating scenarios release SO 2 from sulfides, dominating stratospheric sulfate aerosols and causing long-term global cooling.Conversely, in moderate-temperature heating scenarios, the predominant release of CO 2 from hydrocarbons induces global warming.
In carbonate rocks with small amounts of sulfate, high-temperature heating contributes to global cooling.Locations abundant in evaporites, specifically sulfate rocks, experience more pronounced cooling, surpassing the greenhouse effect of CO 2 .For instance, the end-F LIPs, characterized by abundant sulfate rocks, exhibit significant cooling due to the sunlight-blocking effect of sulfate aerosols.
The asteroid impact at the K-Pg boundary, causing short-term heating, leads to the generation of soot from hydrocarbons in sedimentary rocks.This induces global cooling irrespective of the quantity of sulfate rocks.
The interaction between magmas and sedimentary rocks results in the significant release of SO 2 and CO 2 gases through sill contact metamorphism, which in turn leads to diatreme pipe eruptions 1,34 .Evidence of contact metamorphism has been documented in association with mass extinctions in the end-F, D-C, end-Permian, and end-Triassic, with numerous sills reported across a wide range 1,7,[34][35][36]41,42 (see Table 3). The qantity of sills and sedimentary rocks in contact necessary to trigger global changes is present within each LIP 1,7,[34][35][36]41,42 .Utilizing Arrhenius's equation, the experimental temperature over minutes to one day is converted to a temperature scale spanning 100 years for sills.This temperature correlates with the temperature of the aureole surrounding the sill.
Examining specific mass extinction events, such as the end-Triassic, end-Permian, and Late Devonian, reveals distinct patterns of SO 2 and CO 2 emissions corresponding to low, intermediate, and high-temperature heating, influencing global cooling or warming.The Chicxulub impact at the K-Pg boundary, characterized by shortterm heating, released substantial quantities of SO 3 , soot, and CO 2 gases from sedimentary rocks, contributing to global cooling 15,27 .
Global cooling resulting from a single volcanic eruption or impact event occurs over 2-3 years, mostly recovering within approximately 10 years 15 .Subsequently, a period of warming follows 43 .The prominence of either cooling or warming depends on the heating temperature.Specifically, large-scale volcanic regions with a diameter of approximately 2000 km, containing diverse sedimentary rocks, establish heating temperatures as a key factor in climate change dynamics.The influence of heating temperature on mass extinctions Our analysis underscores the crucial role of sedimentary rock temperature in shaping climate changes during mass extinctions.We have identified distinct temperature-dependent processes that drive global cooling and warming.Low-temperature (400-750 °C) heating of sulfide leads to the release of SO 2 , triggering global cooling through the formation of sulfuric acid aerosols.This phenomenon was observed during sill development in the continental rift volcanism of Central Atlantic Magmatic Province.Central Atlantic Magmatic Province volcanism, which contributed to the end-Triassic mass extinction (Fig. 6).Intermediate temperature (700-1100 °C) heating of hydrocarbons and carbonates releases significant amounts of CO 2 , resulting in global warming.Such occurrences were observed during felsic extrusive volcanic eruptions in the Emeishan and Siberian large igneous provinces, leading to the G-L mass extinction and end-Permian mass extinction, respectively (Fig. 6).Hightemperature (1000-1300 °C) heating produces both SO 2 from sulfate and sulfide, leading to global cooling.Notably, this was observed during ultrabasic kimberlite-like eruptions associated with the onset of the end-Frasnian mass extinction.Furthermore, at medium-high temperatures and short heating durations, the formation of soot from hydrocarbons plays a significant role in global cooling.

Conclusions
This study focused on the mechanisms of climate change during mass extinctions in the Phanerozoic Eon.These variations include cases where global cooling is predominant and cases where global warming is predominant, but the reasons for these are not always clear for each mass extinction event.I showed that the temperature to which sedimentary rocks are heated determines these climate changes.Low-temperature heating of sulfides causes global cooling through the release of sulfur dioxide (SO 2 ), while intermediate-temperature heating of hydrocarbons and carbonates results in the release of large amounts of carbon dioxide (CO 2 ), leading to global warming.Furthermore, high-temperature heating generates SO 2 from sulfates, contributing to global cooling.In the case of large-scale volcanic activity, magma intrudes into sedimentary rocks, generating the aforementioned volatile gases depending on the temperature of contact heating.These gases are then released into the atmosphere through gas explosions, causing the aforementioned climate changes.Additionally, short-term medium to hightemperature shock heating resulting from asteroid impacts generates soot from hydrocarbons, also contributing to global cooling.
To further clarify the mechanisms underlying climate control through the heating temperature of sedimentary rocks during significant mass extinctions, future investigations should prioritize the acquisition of data pertaining to the following: 1. Sea surface temperature (SST) during the end-Devonian mass extinction and end-Jurassic mass extinction (Fig. 6).2. Coronene index and the quantities of five-to seven-ring PAHs during the end-Ordovician mass extinction and end-Jurassic mass extinction (Fig. 5).
By obtaining these additional data sets, we can develop a more comprehensive understanding of the role of heating temperature in driving climate changes during mass extinctions.
This paper has elucidated the causal mechanisms of climate change associated with five well-documented mass extinctions.It has achieved this by clarifying the relationship between heating temperature, the release amount of SO 2 , CO 2 , and soot, which control global climate, and the heating temperature proxy, coronene index, recorded in sedimentary rocks.Additionally, the study reveals the evident connection between climate change itself and its causal locations.

Determination of CO 2 and SO 2 gas release temperature
This study aims to determine the temperatures at which CO 2 and SO 2 gases are released from organic carbon-rich limestone, specifically those containing 1 wt% gypsum (sulfate) and 1 wt% anhydrite (sulfate), respectively.The experiments were conducted using the TG/STA (thermo balance) coupled with a Quadrupole Mass Spectrometer (QMS) system (Fig. 2).The rock sample utilized in the experiments is a Changhsingian limestone sourced from Meishan, China, containing a total organic carbon (TOC) content of 0.5%, supplemented by Miocene gypsum obtained from San Jose, California and artificial anhydrite.The heating temperature during the experiments was gradually increased by 10 °C within a duration of 1 min, up to 1500 °C.This approach allowed for controlled heating of the sedimentary rocks and the measurement and characterization of gas release.The relevant data are presented in Fig. 2. In this study, we also utilize previous gas release measurements, specifically the type and quantity of gas released from heated mudstone and limestone 11 , as well as high-sulfur coal 17 , to determine the temperature at which CO 2 and SO 2 gases are released.

Investigation of soot formation temperature
In addition to the aforementioned experiments, heating experiments were performed on sedimentary rocks to examine soot formation.Raman analyses were employed to identify the presence of soot.Sixteen powder samples were prepared from the previously mentioned limestone (Fig. 3a, b), while eight powder samples were prepared from a Changhsingian black shale from Shangsi, China, with a higher total organic carbon (TOC) content of 2.0% (Fig. 3c).For each rock sample, varying amounts of powder (1, 2, or 5 g) 11 were placed into quartz ampoules, which were then evacuated and sealed at a pressure of 2 × 10 −4 Pa.Each ampoule was individually heated in an electric furnace, with temperatures ranging from 350 to 1100 °C in increments of either 50 °C or www.nature.com/scientificreports/ 100 °C11 .Each maximum temperature is kept for 1 min to 1000 min (Fig. 3a-c).The maximum temperature within ± 10 °C is kept for approximately 5 min in the 1-min case.Five minutes is used as heating time for soot in the 1-min case (Fig. 3).

Arrhenius equation
To determine the frames of experimental data in Fig. 4, we utilized an equation derived from the Arrhenius equation: T 1 = E a T 2 /(E a − T 2 * ln[t 2 /t 1 ]); T 1 : temperature in actual heating time; Ea: activation energy (74 kcal/mol for n-C 16 alkane 44 ; 74 kcal/mol for calcium carbonate 45 [310 kJ/mol = 74 kcal/mol], and 67 kcal/mol for pyrite 46 ; T 2 : temperature in the experiment; t 1 : heating time (s) in the experiment; t 2 : actual heating time 11 .There are no significant differences in amount of hydrocarbon and CO 2 gases generated from under different pressure between 0.1 MPa (1 atm) and 30 MPa (1 km depth underground, which is maximum depth treated in this study) based on experiment data of oil and coal combustion 11,44,47 .Pressure difference on average causes ca. 10 °C higher, which can be neglected 11,44 .

Compilation of SST and coronene index
The analysis in this study focuses exclusively on sea surface temperature (SST) data from low latitudes, obtained from sources [48][49][50][51][52] .It is important to note that surface temperature anomalies at low latitudes consistently exhibit intermediate values, indicating proximity to average values 4 .The error margin for SST anomalies in geological ages is approximately 1 °C, which includes approximately ± 0.5 °C depending on the sample location to obtain the average value and approximately ± 0.5 °C depending on the detection of the largest anomaly for abrupt short-term events from sedimentary rocks, typically deposited at a rate of 1-100 mm kyr −1 , except for impact ejecta sediment 4 .SST anomalies from various geological ages were converted to global surface temperature anomalies and land-surface temperature anomalies using Figure 1d of Kaiho 4 , which was generated from global cooling and warming (recovery) data obtained through climate model calculations 27 .These compiled data are presented in Table 1.
The coronene index, proposed by Kaiho et al. 9 , is defined as the ratio of coronene (seven-ringed polycyclic aromatic hydrocarbons [PAH]) to the sum of benzo[e]pyrene (five-ringed PAH), benzo[ghi]perylene (six-ringed PAH), and coronene itself.It serves as an indicator of heating temperature.Previous studies [9][10][11][12]31 have utilized this index to investigate mass extinction events. Thecoronene index data are compiled in Table 1.
Based on these findings, this paper proposes a theory that seeks to uncover the mechanisms responsible for the observed global cooling and warming during mass extinctions by synthesizing the summarized data.

Figure 1 .
Figure 1.Percentages of species extinction in marine animals and tetrapods during major and minor mass extinctions.Blue columns indicate extinction values for marine species, while red columns represent extinction values for terrestrial species 4 .Open histograms in H-A depict species extinction percentages for 2060-2080 CE 53 .The abbreviations O, F, D, G-L, P, T, J, K-Pg, and H-A correspond to the following time periods: Ordovician, Frasnian, Devonian, Guadalupian-Lopingian transition, Permian, Triassic, Jurassic, Cretaceous-Paleogene boundary, Holocene-Anthropocene.These extinction percentage data are directly comparable due to the use of similar methods, including the conventional method and substage intervals 4 .The numbers 1 to 5 denote the five major mass extinctions.Each silhouette represents a representative vertebrate animal from the respective age.The end-Jurassic mass extinction event is associated with the Morokweng impact crater formation in South Africa 3 .SST: presence of SST data.

Figure 2 .
Figure 2. Temperature-dependent release of CO 2 and SO 2 gases from organic carbon-rich limestone containing 1 wt% sulfate (gypsum and anhydrite).(a) Temperature profiles illustrating the release of CO 2 gas from organic carbon-rich limestone.(b)Temperature profiles illustrating the release of SO 2 gas from samples containing 1 wt% gypsum (sulfate) and 1 wt% artificial anhydrite (sulfate).The curve labeled 44 amu represents the temperature for CO 2 release, while the curve labeled 64 amu (b) corresponds to the temperature for SO 2 release from gypsum and anhydrite.The molar ratio of SO 2 /g gypsum or anhydrite to CO 2 /g limestone is 4%, which corresponds to a weight ratio of 6%.

Figure 3 .
Figure 3. Temperature-dependent release of soot from sedimentary rocks, substantiated by compelling evidence.Temperature profiles illustrate the release of soot from organic carbon-rich limestone (a,b) and organic carbon-rich black shale (c).Refer to the "Methods" section for detailed information on samples and methodologies.During heating in panel (a) at 1050 °C and 1100 °C, the limestone powder underwent partial melting, indicating the conversion of CO 2 to CO and O 2 , accompanied by the decomposition of soot.Importantly, this reaction is improbable in an open system.An examination of the Raman spectrum of a black material found on the inner wall of the ampoule confirmed its composition as soot (d).The outer diameter of the ampoules in panels (a,b) and panel (c) is 25 mm and 30 mm, respectively.

3 Figure 4 .
Figure 4.The temperature and heating durations required for the production of SO 2 , CO 2 , and soot were determined by integrating data from Figs.2 and 3, along with findings from published papers.The figure illustrates the formation processes of various gases: SO 2 from sulfide (based on Wang et al.17 using highsulfur coal, with pyrolysis experiments conducted for 15 min at a rate of approximately 30 °C per minute, and Kaiho et al.11 using a heating rate of 8 °C per minute up to 1100 °C), CO 2 from hydrocarbon and carbonate (derived from Kaiho et al.11 and the present study), SO 2 from sulfate (from the present study), and soot from hydrocarbon (based on Krestinin et al.19 [isothermal pyrolysis of acetylene, benzene, and diacetylene], Yoshihara et al.20 under a shock-heated condition, He et al.21 under 2 and 4 atm shock temperature and pressure, model calculation of Chen et al.22 , and the present study).The data points are represented by red and right-side blue bars, and black dots for 1-min heating, all obtained from the present study.The term "Low oxygen" corresponds to a pressure of 2 × 10 -4 Pa.The colored areas on both sides of the frame, which encompass the data plots, represent the fitting of the experimental data to the Arrhenius equation.The Arrhenius equation elucidates the relationship between production temperature and heating time.The presence of experimental data points within these colored areas indicates how well they conform to the predicted behavior based on the Arrhenius equation.The production temperature on the horizontal axis indicates the heating temperature in each experiment.

Figure 5 .
Figure 5.The three classifications of mass extinctions depicted in the upper part of this figure are based on surface temperature anomaly and the Coronene index (a sedimentary rock heating temperature index).The lower part of the figure illustrates the heating temperature based on experimental data and the major emission gases/particles influencing climate change.The colored vertical bars in the upper figure indicate the average values of the Coronene index for each mass extinction.Abbreviations and colors represent different events: End-F (end-Frasnian): green dots, G-L (Guadalupian-Lopingian transition): orange dots, End-P (end-Permian): red dots, End-T (end-Triassic): blue dots, and K-Pg (Cretaceous-Paleogene boundary): black dots.Specific emission gases in the lower figure, such as SO 2 , CO 2 , or soot, are determined by the heating temperature.The relationship between the three temperature scales connecting the upper and lower figures is based on the Arrhenius equation11 .The upper temperature scale indicates the instantaneous heating temperature, such as that resulting from impact and experiments on the formation of PAHs, including coronene.While the temperature at the impact point exceeds 10,000 °C, the heating temperature decreases to a few hundred degrees moving from the center of the crater outward.The middle temperature scale represents the temperature during 5 min of heating based on Kaiho et al.11 and experimental data from this paper.The bottom temperature scale represents the temperature heated by sill according to Aarnes et al.18 , lasting for 100 years.

Table 3 .
Lithologies of large igneous provinces for five mass extinctions and asteroid impact target rocks.